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The expressions for average densities of currents and charges induced by a weak electromagnetic field in 
spatially inhomogeneous systems are obtained. The case of finite temperatures is considered. It is shown 
that average values are separated into "basic" and "additional" parts. The former depends on electric fields, 
and the latter depends on derivatives of electric fields on coordinates. Semiconductor quantum wells, wires 
or dots may be considered as spatially inhomogeneous systems. 



l.R. J. Kubo [1] obtained the formula for the conduc- 
tivity tensor CTq,/3(w), applicable in a case of spatially 
homogeneous systems and electric fields E(t) indepen- 
dent of spatial coordinates. This formula takes into ac- 
count exactly the interaction of current carriers with the 
medium. Consequently, it is a powerful tool for the solu- 
tion of concrete conductivity problems in solids. 

Currently there is a great interest in experimental and 
theoretical study of low-dimensional semiconductor ob- 
jects such as quantum wells, wires and dots. Because of 
this, the generalization of the Kubo formula in a case of 
spatially inhomogeneous systems and spatially inhomo- 
geneous fields is urgent. Kubo [1] used the operator of 
interaction of current carriers with the electric field as 



Uk = - ^e,r,E(t), 



(1) 



where ei and r.^ are the charge and radius - vector of the 
i -th particle, E(t) is time-dependent, but homogeneous 
in the space electric field. However, it follows from the 
Maxwell equations, that the time-dependent electric field 
necessarily depends on coordinates, so the use of Eq. (1) 
is always a certain approximation, if E depends on t. 
Our task consists of obtaining the additional terms in the 
Kubo formula, which contain derivatives of the electric 
field on coordinates, as well as taking into account the 
heterogeneity of physical systems. 

2. Let us consider a system of N particles of charge 
e and mass m in a weak electromagnetic field, charac- 
terized by intensities E(r,t) and H(r, i). We introduce 
vector A(r, and scalar ip{'R.,t) potentials in terms of 
which the fields are expressed (the gage of potentials is 
arbitrary) . For the sake of generalizing our task, we shall 
suppose that the system of particles is placed in a con- 
stant magnetic field H. The vector potential -4(r) corre- 
sponds to H. The total Hamiltonian Htotai is as follows 



where 



n 



^ 2n 



+ y(ri...rAf), 



(3) 



Pi ^P,~ {e/c)A{ri), Pi = -ih[d/dvi), V(ri...rAr) is 
the potential energy, including interaction between par- 
ticles and external potential. U is the interaction energy 
of particles with the electromagnetic field 



U = Ui + C/2, 



Ui^~^J dVj(r)A(r,i) + J (frp{r)ip{Y,t), 



2mc 



(4) 



(5) 



(6) 



n 



total 



(2) 



In Eqs. (5), (6) the operators of current density 



2 m 
and charge density 

= ^Pi(r), pi{r) = e^(r - r,) 

i 

are introduced. The operator U2 is out of the framework 
of the linear approximation on the electromagnetic field. 

It is supposed that on the indefinitely removed dis- 
tances there are no charges and currents, and also that 
on times t — + —00, the fields E(r, t) and H(r, t) arc equal 
to 0, which corresponds to adiabatical switching on the 
fields. Generally, the operator U of the interaction of 
particles with the electromagnetic field is not expressed 



1 



through fields E(r,t) and H(R, t). Accordingly, the op- 
erators ji{r,t) and pi{r,t) are expressed only through 
potentials A(r, t) and <p(r, t) as follows 

e I /■* 

jia(r,i) = p{r,t)A^{r,t) + - / dt'[Ui{t'),j^{r,t)], 

TnC II I _ryr^ 



(Pi(r,i))aB/ar = 



Tie 



dt' 



Pi{v,t)='-j dt'[Ui{t'),p{v,t)], 



(7) 



where subscript 1 means linear approximation on the 
fields. [F, Q] = FQ — QF is the commutator of oper- 
ators F and Q. 

However, all the observable values must be expressed 
through fields. This also is true for average values 
(iia(r, t)) and {pi{r,t)), where the averaging is intro- 
duced by 



<...>= 



Sp{e^p{-I3n) ...} 
Sp{eM-Pn)} ' 



kT' 



(8) 



The problem lies in what kind of interaction to use a 
containing vector potential or an electric field, which was 
discussed earlier in [2] with reference to a task of the light 
dispersion in bulk crystals. However, since we must solve 
other problems with spatially inhomogeneous systems, it 
is necessary to come back again to this topic. 

3. For the case T = average values of densities of cur- 
rent and charge, induced by external fields in a spatially 
inhomogeneous medium, are calculated in [3]. The ob- 
tained expressions contain the symbol (0|...|0) , where |0) 
is the wave function of the ground state of the system. 
By replacing averaging (0|...|0) with averaging Eq. (8) 
[3] we obtain 



(ilc<(r,i)> = {hoc{r,t)) E + {jla{T^,t)) OE/dr 



(9) 



[pi{r,t)) = {pi{r,t)) E + {piir,t)) sE/ar, (10) 



where 



(ii„(r,t))B = J J d\' j' dt' 
x([j„(r,t), d0{v',t')])E0{v',t'), 

{pi{r,t))E = ^J d^'r' j' dt' 
x{[p(v,t), dp{v' ,t')]) Ep{v' ,t'), 



' mc or a 



(11) 



(12) 



I 

he 



J d^r' J dt'{[Ur,t),Y(,^{r',t')]) 
^da/3{r',t') 



(13) 



x([p(r,t),y,.(r',0])^^^^, (14) 



d(r)=rp(r), ^^^(r) = r^i^(r), (15) 



a(r,t) = -c I dt'E{r,t'). (16) 

The correctness of the replacement of averaging 
(0|...|0) with ( ...) proves to be true in the further com- 
parison of our results with the results of [4] and [5]. 

4. We transform Eqs. (9) - (14) so that at transition to 
a limit E(r, t) ~ E(t), the Kubo formula for conductivity 
[1] will be obtained. We use a ratio [1,4] 



h 



{{m,Q{t')]) 



dX{^^F{t + zhX)), (17) 



true for any pair of operators F and Q, and wc also use 
the integration on variable t' and r' in parts. It results 
in 



(.7i„(r,i)) =(ji„(r,i))(i) + (.n„(r,0) 



(2) 



(18) 



where the first part contains an electric field, and the 
second contains its derivative on coordinates, i.e. 



(ii„(r,t))W= / /* dt' f d\ 

J J—oo Jo 

X { Mr', t') jcc (r, t + inX)) Ef) (r', t'), 



(19) 



(i.„(M))(^) = ^(.,(r))^4^ 
mc or alpha 

~ J d'r' J^' dX{Y0^{r')j^{r,ihX)) 

^ dapir'.t) 
dr'^ • 

Similarly, we obtain 

(pi(r,i)) =(pi(r,i))W + (pi(r,i))(2), 

(/)i(r,t))(i) = / d^r' f dt' f dX 

J J-oo Jo 

x{jpir',t')pir,t + ihX))E0{r',t'), 



(pi(r,t))(2)=_i / d'r' r dX 



(20) 



(21) 



(22) 



x(y,,(r')p(r,aA))^^^ 



(23) 
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Let us emphasize that the separations of the average val- 
ues Eqs. (18) and (21) in two parts does not coincide 
with the separations Eqs. (9) and (10), convenient only 
at T = 0. 

The contributions with the superscript (1) contain an 
electric field. We shall name them "basic". The contribu- 
tions with the superscript (2) contain derivatives of the 
electric field on coordinates. We shall name them "addi- 
tional" . The reason for introducing basic and additional 
contributions is that for the solution of any problem, in 
which the field E(r, t) is spatially inhomogeneous, it is 
possible to estimate magnitudes of the additional contri- 
butions and to determine whether it is necessary to take 
them into account or to reject them. 

5. We compare our results with the conclusions of [5], 
where (jia(r,t)} is also calculated at finite temperatures. 
This expression differs from ours which only contains an 
electric field. Generally speaking, integrating on r' in 
parts, it is possible to get rid of derivative dEp(r' , t)/dri^, 
proceeding to the formulas which only contain fields. 
However, the opposite procedure, transition from fields 
to derivatives, is not always possible in basic contribu- 
tions where fields are always kept. 

The separation of Eq. (18) in two parts coincides with 
the separation of [5]. The expressions for the basic con- 
tribution with superscript (1) coincide. But instead of 
Eq. (20) in [5] for the additional contribution, the result 



dt 



is obtained, where 



= - j dV ^„^(r,r',0) 
x[E0{v' ,t) - Ep{v,t,)\ (24) 



¥5a/3(r,r',T) = / d\{ji3{r')j^{v,T + ih\)). (25) 
Jo 

Integrating both parts of Eq. (24) on time, we have 
(ji«(r,t))(2)=--£_(p(r))a«(r,t) 

dh' j\\{30{v')3cc{v,in\))ap{v',t). (26) 
Similarly, we obtain 

(p(r,i))(2) = l| dS' j\x 

{jp{v')p{v,ih\))afi{v\t). (27) 

We have verified that it is possible to proceed from our 
results in Eqs. (20) and (23) to the results of Eqs. (26) 
and (27), integrating on r' in parts. 

6. Eqs. (20) and (23), unlike Eqs. (26) and (27), 
contain operators rj of coordinates of particles. Indeed, 
the definitions in Eqs. (15) may be written as 



d(r) = e^ Ti5{Y-ri 



(28) 



(29) 



But the average values (ji(r,t)) and (/9i(r,t)) should 
not depend on the point of readout of the coordinates rj. 
This means that Eqs. (20) and (23) contain only non- 
diagonal elements of the operators r,, and the diagonal 
elements can be excluded. 

To demonstrate this statement let us transform Eq. 
(29) so that the operator stands only on the left 

^/37 (r) = - 1^ hi P(r) + X! ^^1^ ^i-r ('^) • (^0) 



Substituting Eqs. (28) and (30) in Eq. (20), we obtain 



+ 



mc 



2mc 



J d\' ^ dx{p{v')j„{v,inx)) 

xdiv a(r',t) 

[ d^r' [ dXy2{ril3jij{r')ja{r,ihX)) 
cj Jo . 

aa^(r^^) 

Let us separate the operator in two parts: 



(31) 



(32) 



where superscripts d and nd mean diagonal and non- 
diagonal contributions, respectively, and 



{n\rf\m) = (n|ri|m)(n|m). 



(33) 



We then substitute Eq. (32) in the last term of Eq. 
(31) from item 4. Taking advantage of commutativity 
[H, rf] = 0, it is possible to show that 

"c / d'r' jyxY,{rt0jn{r')ja{r,inX)) 



dai3{r',t) e ^ 
^;:^ = -;;^E(-^/3<5(r-r,)) 



dr„ 



\ j d:'r' jyxY,{ki{r:')U^,inX)r^^{inX)) 

da0{v',t) 



(34) 



Using Eqs. (32) and (34), we obtain from Eq. (31) ex- 
pressions, which do not containing the operators rf, 
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(iia(r,t)) 



(2) 



2mc , 

x{p{r')ja{r,ihX)) divsiivtj) 



jUr')j^{r,inX) 



dai3{r',t) 
dr', 



(35) 



Similarly, for the additional contribution in average in- 
duced density of charge we have 



ih 



(P(r,t))(^) = ^ / rfV / dX 
X ( p(r') p(r, ihX)) div a(r', t) 

dxY,{3^',{Ap{^,inx)r2^(inx) 



^,%,ir')pir,inX))^^^. 



(36) 



7. We make the Fourier-transformation of the electric 
field as follows 

E^{r,t) = Ei+\r,t) + Ei-\r,t), 
Ei+\r,t) = ^ J d^kj^ du;£;„(k,a;)e^('^'— *), 

4-)(r,i) = (i?(+)(r,t))*, 



d^r / dtE^{v,t)e-'^^''-'^'\ 

J —OO 



(37) 



The average induced density of current may now be writ- 
ten as 



(iia(r,i)) 



(2 



J Jo 



djjj 



xaa/3(k,a;|r)£;3(k,w)e'('''-'"*) + ex., 



(38) 



where crQ/3(k, w | r) is the conductivity tensor, dependent 
on spatial coordinates. Using Eqs. (19) and (35) for 
basic and additional contributions to the average induced 
density of current, we obtain 



(7„/3(k, w I r) = (T^j (k, w I r) -I- a^j (k, w | r), 
where 



(39) 



d'r' 



dt 



x{jp{v-v',-inX)Uv,t))e-'^'^^'- 





dX 

iiot) 
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(40) 



/3 



dX 



^ / d'r'e-'^^' 
2muj J 

x{p{r-r',-inX)ja{r)) 



As an approximation, when the electric field does not 
depend on spatial coordinates, i.e. E(r, t) ~ E(t), we 
make the Fourier-transformation 

1 f°° 

{jia{r,t))h ^ 2^ ^ 



xaaf3{i.o\r) Ef3 (lo) e~ 



+ C.C.. 



(42) 



where subscript h means a field, homogeneous in space. 
It is easy to see that 



aapiuJ I r) = (Jo3 (k = 0, w I r). 
Then with the help of Eqs. (40) - (42) we obtain 
r0 



(43) 



0-a/3('^|r) = 



dt 



dX{j0{-inX)Uv,t))e'^\ (44) 



where the current operator is 
■4 = e ^ : 



(45) 



Eq. (44) is the generalization of the Kubo formula in 
the case of a spatially inhomogeneous medium, where 
the conductivity tensor depends on r. In the case of a 
spatially homogeneous medium, the tensor aajj{ui\v) does 
not depend on r and we obtain from Eq. (44) 



^aM'^) =vj^ '^^ X dX{j0{-ihX)J^{t))e'^\ 



(46) 



^ J Jo ^ 

-r^\-ihX)ji^{r - r', -i?iA)j„(r)). (41) 



where Vq is the normalization volume. The obtained for- 
mula coincides with the Kubo result [1] if we substitute 
oj for —oj on the RHS and take into account that in [1] 
the volume Vq = 1. 

It follows from Eq. (42) that the additional contribu- 
tions to the average density of the induced current re- 
sult in additional contributions to the conductivity ten- 
sor which contains the factor k-.,/uj. If a field E(r,t) is 
the plane wave (at a monochromatic irradiation) or wave 
package ( at a pulse irradiation). A; ~ w/c and additional 
contributions contain a small factor v/cin comparison to 
the basic contributions, where v is the particle velocity 
in the system. 

However, this estimation is not always correct for spa- 
tially inhomogeneous systems (for example for semicon- 
ductor quantum wells, wires or dots). It is possible to 
consider the field E(r, f) as external, or as a stimulating 
field, only when we calculate the densities of the induced 
current and charge in the lowest order on interaction of 
the field with the system of charged particles. Such an 
approximation is acceptable in the case of a quantum well 
under the condition [6,7] 

7r < 7, 
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where 7, (7) is the radiative (non-radiative) broadening 
of electronic excitation. 

Otherwise, when 7^ 7, it is necessary to take into 
account the interaction of the field with particles in all 
orders of a perturbation theory. Then the field E(r, t) 
is the genuine field within the limits of low-dimensional 
objects. 

This field cannot yet be represented as a superposi- 
tion of plane waves, for which k = u)/c. For example, 
the genuine field strongly varies inside of a quantum well 
along the z axis which is perpendicular to the plane of 
a quantum well. The light falls along the z axis and 
the frequency uj is in resonance with one of the excita- 
tion energy levels of the electron system in a quantum 
well belonging to the discrete spectrum [8,9]. Then the 
values kd ~ 1 are essential, where d is the width of a 
quantum well. Instead of a small factor v/c, we obtain a 
factor 



If the wave length is A ^ d, it may appear that the new 
factor M is much greater than v/c. In concrete cases, it is 
necessary to estimate its magnitude. In [8,9] we assumed 
that M <C 1 and neglected the additional contributions 
to average values of the induced densities of currents and 
charges (it was considered in the case T = 0). 

Thus, Eqs. (40) - (42) represent the analogue of the 
Kubo formula in a case of spatially inhomogeneous elec- 
tric fields and systems. The average values of the induced 
charge density arc also calculated in this case (see Eqs. 
(21), (22) and (37)). 

The problem of the transition from initial expressions 
(containing vector and scalar potentials) to expressions 
for average values (containing electric fields and deriva- 
tives from field on coordinates) is solved. 

Substituting the obtained expressions for the average 
densities of current and charge in the Maxwell equations, 
it is possible in principle to determine the genuine fields 
inside and outside of semiconductor low-dimensional ob- 
jects. Thus, it is possible to calculate the light reflection 
and absorption coeflacients of these objects (see, for ex- 
ample, [8,9]). 
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